FFT_2D Function

public function FFT_2D(signal, method, threads) result(result)

Perform a 2D Fast Fourier Transform on a signal

This function takes a signal and performs a Fourier Transform using the specified method. The available methods are:

  • "NAFPack_FFT_2D": Fast Fourier Transform using NAFPack
  • "FFTW_FFT_2D": Fast Fourier Transform using FFTW
  • "FFTW_FFT_2D" + threads: Fast Fourier Transform using FFTW with multithreading

Arguments

Type IntentOptional Attributes Name
complex(kind=dp), intent(inout), DIMENSION(:, :) :: signal
character(len=*), intent(in) :: method
integer, intent(in), optional :: threads

Return Value complex(kind=dp), DIMENSION(SIZE(signal, 1) ,SIZE(signal, 2))


Calls

proc~~fft_2d~~CallsGraph proc~fft_2d FFT_2D fftw_cleanup_threads fftw_cleanup_threads proc~fft_2d->fftw_cleanup_threads fftw_destroy_plan fftw_destroy_plan proc~fft_2d->fftw_destroy_plan fftw_execute_dft fftw_execute_dft proc~fft_2d->fftw_execute_dft fftw_init_threads fftw_init_threads proc~fft_2d->fftw_init_threads fftw_plan_dft_2d fftw_plan_dft_2d proc~fft_2d->fftw_plan_dft_2d fftw_plan_with_nthreads fftw_plan_with_nthreads proc~fft_2d->fftw_plan_with_nthreads

Source Code

    FUNCTION FFT_2D(signal, method, threads) RESULT(result)

        COMPLEX(dp), DIMENSION(:, :), INTENT(INOUT) :: signal
        CHARACTER(*), INTENT(IN) :: method
        INTEGER, OPTIONAL, INTENT(IN) :: threads
        COMPLEX(dp), DIMENSION(SIZE(signal, 1) ,SIZE(signal, 2)) :: result

        IF(method == "NAFPack_FFT_2D")THEN
            result = NAFPack_FFT_2D(signal)
        ELSE IF(method == "FFTW_FFT_2D" .AND. .NOT. PRESENT(threads))THEN
            result = FFTW_FFT_2D(signal)
        ELSE IF (method == "FFTW_FFT_2D" .AND. PRESENT(threads))THEN
            result = FFTW_FFT_2D_threads(signal, threads)
        ELSE
            STOP "ERROR : Wrong method for FFT_2D"
        END IF

    END FUNCTION FFT_2D